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C — , I We consider the open spin-s XXZ quantum spin chain with N sites and gen- 

eral integrable boundary terms for generic values of the bulk anisotropy parame- 

. | ter, and for values of the boundary parameters which satisfy a certain constraint. 

We derive two sets of Bethe Ansatz equations, and find numerical evidence that 
together they give the complete set of (2s + 1)^ eigenvalues of the transfer ma- 
trix. For the case s = 1, we explicitly determine the Hamiltonian, and find an 
expression for its eigenvalues in terms of Bethe roots. 
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1 Introduction 



Much of the impetus for the development of the so-called Quantum Inverse Scattering 
Method and Quantum Groups came from the effort to formulate and solve higher- 
spin generalizations of the spin-1/2 XXZ (anisotropic Heisenberg) quantum spin chain 
with periodic boundary conditions p]-[5]. (See also [6] and references therein.) For 
the open spin-1/2 XXZ chain [7], Sklyanin [5] discovered a commuting transfer matrix 
based on solutions of the boundary Yang-Baxter equation (BYBE) [9]. This made it 
possible to generalize some of the above works to open spin chains [ID]- [IE]- Although 
general (non-diagonal) solutions of the BYBE were eventually found [T7J [THJ, [15] , Bethe 
Ansatz solutions were known only for open spin chains with diagonal boundary terms. 
Further progress was made in [191 [20], where a solution was found for the open spin-1/2 
XXZ chain with general integrable (i.e., including also non-diagonal) boundary terms, 
provided the boundary parameters satisfy a certain constraint. It was subsequently 
realized [21] that two sets of Bethe Ansatz equations (BAEs) were generally needed in 
order to obtain the complete set of eigenvalues. (Recently, the second set of BAEs was 
derived within the generalized algebraic Bethe Ansatz approach of [19] by constructing 
[22] a suitable second reference state. For other related work, see e.g. [23]- [31], and 
references therein.) 

We present here a Bethe Ansatz solution of the iV-site open spin-s XXZ quantum 
spin chain with general integrable boundary terms, following an approach which was 
developed for the spin-1/2 XXZ chain in [M], and then generalized to the spin-1/2 XYZ 
chain in [35J. Doikou solved a special case of this model (among others) in [32] using 
the method in [20], and also computed some of its thermodynamical properties. She 
further investigated this model in [33] using the method in [19] . We find a generalization 
of Doikou's constraint on the boundary parameters, as well as a second set of BAEs, 
which we argue is necessary to obtain the complete set of (2s + 1)^ eigenvalues. 

The outline of this paper is as follows. In Sec. [2] we review the construction of the 
model's transfer matrices, and the so-called fusion hierarchy which they satisfy. In Sec. 
[HI we identify (as in |34j) the Q operator with a transfer matrix whose auxiliary-space 
spin tends to infinity, and obtain a T — Q equation for the model's fundamental (spin- 
1/2 auxiliary space) transfer matrix. Using some additional properties of this transfer 
matrix, we arrive at a pair of expressions for its eigenvalues in terms of roots of cor- 
responding BAEs, as well as a corresponding constraint on the boundary parameters. 
We provide numerical evidence of completeness of this solution for small values of s 
and N. In Sec. H]we explicitly determine the Hamiltonian for the case s = 1, and find 
an expression for its eigenvalues in terms of Bethe roots. We end in Sec. [5] with a brief 
further discussion of our results. 
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2 Transfer matrices and fusion hierarchy 



Sklyanin [8] showed for an iV-site open XXZ spin chain how to construct a commuting 
transfer matrix, which here we shall denote by t^'^\u), whose auxiliary space as well 
as each of its N quantum spaces are spin- 1/2 (i.e., two-dimensional). In a similar way, 
one can construct a transfer matrix t^' s \u) whose auxiliary space is spin-j ((2j + 1)- 
dimensional) and each of its N quantum spaces are spin-s ((2s + l)-dimensional), 
for any j,s G {|,1, §,...}. The basic building blocks are so-called fused R and K T 
matrices. The former are given by [36], [2]- [5] l. 



{a}{b}\ 



2j 2s 



) = P {a} P { b } II II } ^ + j — S — 1)77) P{ a }P\ 



{b} 



(2.1) 



k=l 1=1 



where {a} = {a±, . . . ,a 2 j}, {b} = {b\, . . . ,62s}, and R^'^\u) is given by 

\ 



( sh(w + 77) 

shu sh?7 

shr/ shu 

\ sh(u + r/) / 



(2.2) 



where rj is the bulk anisotropy parameter. The R matrices in the product ( 12.11) are 
ordered in the order of increasing k and I. Moreover, Pi is the symmetric projector 



P 



(2.3) 



k=i \i=i 



where V is the permutation operator, with V ah ,a k 
example, for the simple case (J, s) = (1, |), 



1; and similarly for PXy- For 



R 



(1.3) 



u 



Pa a RTh'iu - -77) R { M\u + -77) R 
a i, a 2 0.1,0 V O a2,t> V 2 



The fused R matrices satisfy the Yang-Baxter equations 



+ 

ai,<J2 



R 



{u)R 



(M , 
WW 1 



R { t s L(v)R U ' s) 



b {b}{cV 



{a}{c} 



u)R 



{a}{b} 



(2.4) 



(u-v). (2.5) 



We note here for later convenience that the fundamental R matrix satisfies the unitarity 
relation 



'^(—u) — —^(u)l, £(u) = sh(u + 7]) sh(u — 77) . (2.6) 

1 Our definitions of fused R and K matrices differ from those in [33] by certain shifts of the 
arguments. 
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The fused K matrices are given by [TOj [TTJ [TJ] 



p u n 

k=l 



k-1 



lRi] a l\2u+(k + l-2j-l) V ) 



1=1 



x K;fi ) (u+(k-j-l)rj)}P+ } 



(2.7) 



where K ^\u) is the matrix [T71 [18] 
2 (sh a_ ch /?_ ch it + ch a_ sh /?_ sh u) 



sh2w 



e 9 ~ sh 2u 

2 (sh a_ ch /?_ ch m — ch a_ sh /3_ sh «) 



(2.8) 



where a_, /3_ and are arbitrary boundary parameters. The products of braces {. . .} 
in (12.71) are ordered in the order of increasing k. For example, for the case j = 1, 

(A), 1 , _ri.i<>,_ . ^-d), 1 



0,1,0.2 V 

The fused A - matrices satisfy the boundary Yang-Baxter equations [9] 



The fused A" + matrices are given by 



'{a} I 



/<%) 

where the normalization factor 



A" 



(a_,/9_,0_)->(-a+ -/9 + ,e+) 



2j-l I 

f {3 \u) = n ri[-^ 2n + c + k + 1 - 2 ^»] 

1=1 k=l 



is chosen to simplify the form of the fusion hierarchy, see below (12.17p . 
The transfer matrix t^' s '(u) is given by 



u 



trw K^'(u) T&» A"-( M ) T { ^(«) 



where the monodromy matrices are given by products of N fused R matrices, 



>0» i 

l {a},{bl N \} [ 



>(?» , 
R {a},{bW}( U ) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



These transfer matrices commute for different values of spectral parameter for any 
j,f G {-, 1, §, . . .} and any s G {§, 1, |, . . .}, 



t 0> s \u) ,t( j '> s \u') 



0. 



(2.15) 



These transfer matrices also obey the fusion hierarchy [TOl [TTJ [H] § 

t (j-^\u-jr,)t& s \u) = t^ s \u - (j - |)77) + 8^{u)t^ s \u - (j + ifo) , (2.16) 

j = 1, |, . . ., where t^ '^ = 1, and 5^(u) is a product of various quantum determinants, 
and is given by 



u 



2s-l 1 

[ £(u + (a - * - -)77) 



2 

k=0 



2N 

sh(2w - 2r/) sh(2w + 2rj) 



sh(2w — rj) sh(2w + 77) 
x sh(u + q_) sh(u — a_) ch(u + ch(u — (3J) 
x sh(w + a+) sh(u — a+) ch(u + f3 + ) ch(w — . (2-17) 

We remark that the normalization factor f^\u) (12.1 2p has been chosen so that the 
LHS of (12~T6|) has coefficient 1. 

In the remainder of this section we list some important further properties of the 
"fundamental" transfer matrix t^' s \u), which we shall subsequently use to help deter- 
mine its eigenvalues. However, it is more convenient to work with a rescaled transfer 
matrix t^' s \u) defined by 

tt*\u)= (1 * (2-18) 

where 

2s-l . 

g^' s) (u)=Y[sh{u + (s-k + -) V ) (2.19) 

k=l 

(which has the crossing symmetry g^' s \—u — r]) = ±5^2 ^ (u)) is an overall scalar factor 
of the fused R matrix R^ ,s \u). In particular, the rescaled transfer matrix does not 
vanish at u = when s is a half-odd integer. 

This transfer matrix has the following properties: 

t {1 2' s \u + in) = t {1 ?' s) (u) {in - periodicity) (2.20) 

$h s )(- u -r])= t^' s \u) (crossing) (2.21) 



t^> s \0) = -2 3 sh 2Ar ((s + -)rj) chr/sha_ ch/?_ sha+ch/3+1 (initial condition) (2.22) 

2 



2 The derivation of this hierarchy relies on some relations which, to our knowledge, have not been 
proved. See the Appendix for further details. 
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77=0 



= 2 3 sh 2N u — sh «_ ch /3_ sh a + ch /3 + ch 2 u 
+ cha_ sh/3_ ch a + sh /3 + sh 2 « 

- ch(#_ -#+)sh 2 -uch 2 it 1 (semi-classical limit) (2.23) 



22N+1 



D {2N+A)u+{N+2)r) 



ch(6L — 6 + )l for m — > +00 
(asymptotic behavior) 



(2.24) 



As is well known, due to the commutativity property (I2.15p . the corresponding si- 
multaneous eigenvectors are independent of the spectral parameter. Hence, the above 
properties (12.201) - f)2.24p hold also for the corresponding eigenvalues A^' s \u). 



3 Eigenvalues and Bethe Ansatz equations 

We now proceed to determine the eigenvalues of the fundamental transfer matrix. 
Following [3lj, we assume that the limit 

Q{u) = lim t {j -^ s) (u - jrj) (3.1) 
j"-+oo 

exists. § We then immediately obtain from the fusion hierarchy (12.161) an equation of 
the T — Q form for the fundamental transfer matrix t^' s \u), 

Q{u) t^ s \u) = Q(u + V ) + 5 {s \u) Q(u - rj) . (3.2) 

We further assume that the eigenvalues of Q(u) (which we denote by the same symbol) 
have the decomposition Q(u) = f(u) Q(u) with 

M 

Q(u) = Y\ sh(w — Vj) sh(w + Vj + 7]) , (3.3) 
i=i 

which has the crossing symmetry Q(—u — 77) = Q(u). Here M is some nonnegative 
integer. It follows that the eigenvalues A^> s )(u) of t^' s ^(u) are given by 

A^ {u) = Hl{u )9t±I!l + H 2{ u)9^1, (3.4) 
Q{u) Q(u) 

where H\{u) = f(u + 7])/ f(u) and H 2 (u) = S^ s \u)f(u — 77)/ f(u); and therefore, 

H 1 (u-r ] )H 2 (u) = 5 {s) (u). (3.5) 



3 For a more rigorous and extensive discussion of Q operators for the six- and eight-vertex models 
with periodic boundary conditions, as well numerous earlier references, see [37] . 
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The crossing symmetry (I2.2ip together with (13. 4p imply that 

H 2 (u)=H 1 (-u-r ] ). 

We conclude that Hi(u) must satisfy 

- v) Hi(-u -7]) = 5 {s) (u) , 

where 5^(u) is given by (12.171) . 

A set of solutions of (13. 7p for Hi{u) is given by 



(3.6) 



(3.7) 



H ( L ± \u\e u e 2 ,e 3 ) = -2 2 e 2 



2s-l . 

[ sh(u + (s - k - -)ri) 



k=0 



2N 



sh(2w) 



sh(2w + 77) 



x sh(w ± a_ + rj) ch(w ± ei/3_ + 77) sh(-u ± e 2 a + + 77) ch(ii ± e 3 /3 + + 77) , (3.8) 

where ei,e 2 ,e 3 can independently take the values ±1. It follows from (13.61) that the 
corresponding H 2 (u) functions are given by 



H^\u\e u e 2 ,e 3 ) = -2 2 e 2 



2s-l ^ 

[ sh(u + (s - k + -)?]) 



n 2V 



fc=0 



sh(2u + 277) 
sh(2u + 77) 



x sh(w =1= a_) ch(w =f ei/5_) sh(w =f e 2 a + ) ch(w =f €3/?-) 



(3.9) 



An argument from [M] (which makes use of the periodicity (I2.20p ) can again be used 
to conclude that the eigenvalues can be uniquely expressed as 



A^(u) = aH{ ; (n|ei,e 2 ,e 3 )— — h aH^ '{u\ei, e 2 , e 3 )- 



Q(u) 



Q{u) 



(3.10) 



up to an overall sign a = ±1. Noting that the functions Hi ( 13. 81) and H 2 (13.90 have 
the factor ^< s )(w) 2iV (J2HSD in common, we conclude that the eigenvalues of i^' s ">(u) 
(12. 181) are given by 



A ( 2' s) (m) = >|ei, e 2 , e 3 ) v „ , " + ^ ("Id, e 2 , e 3 )- 



Q(u) 



Q(u) 



(3.11) 



where 



#1 («|ei, e 2 , e 3 ) 



#2 (w|ei»e2,e 3 ) 



-2 2 e 2 sh 2 > -(s- hr,) ^ sh(u ±a„ + r 1 ) 
2 sh{2u + 77) 

x ch(u ± ei/5_ + 77) sh(u ± e 2 a + + 77) ch(u ± e 3 /3 + + 77) , (3.12) 
= -2 W « + (s + i),)||^sh( MTa _) 



x ch(u =f ei/3_) sh(n =f e 2 a + ) ch(u =f e 3 /5 + ) . (3.13) 
We have fixed the overall sign a = +1 in (13.111) using the initial condition (I2.22p . 
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The expression (l3~TTj) for A(2> S )( M ) is also consistent with the asymptotic behavior 
(12.2411 if the boundary parameters satisfy the constraint 

a_ + e 1 /3_ + e 2 a + + e 3/ 3 + = eo(0_-0 + )+77A; + -(l-e2)i7r mod(227r) , (3.14) 
where also e = ±1; and if M (appearing in the expression (13.31) for Q{u)) is given by 

M = sN (3.15) 

The requirement that M be an integer evidently implies that 

1 k 

sN — - =)= — = integer . (3.16) 

In particular, for s an integer, k is an odd integer; and for s a half-odd integer, k is 
odd (even) if N is even (odd), respectively. 

Finally, the expression fl3~TTD for A^> s \u) is also consistent with the semi-classical 
limit (12.231) if the {ej} satisfy the constraint 

ei e 2 e 3 = +l. (3.17) 

To summarize: if the boundary parameters {a±, /3±,6±) satisfy the constraints 
(13. 14jl . (13.161) and (13. 1 T[l for some choice (±1) of {e^} and for some appropriate value 
of k, then the eigenvalues of i^' s \u) (I2.18P are given by 

i(J^ ( „ ) = ^ (B | ei ^ | ^ ) 5g^ + flW (B | ei|6l ^ ) 3^!rt, (3.18) 

where H^\u\ei, e 2 , e 3 ) and H^\u\e\, e 2 , e 3 ) are given by (13. 12j) . (13.131) . and 

M<±) 1 k 

Q(±)(n) = ] [ sh(n - vf>) sh(u + vf ] + rf) , M (± ' = S N - - t - ■ (3.19) 

The parameters {w^} are roots of the corresponding Bethe Ansatz equations, 

H { ±\vf ) \e u e 2 ,e z ) Q (±) K (±) +r/) 

' ( L = — > m -, i = l,...,MW; (3.20) 

iyf ) H ±) -7 7 |e 1 ,e 2) 6 3 ) Q^(vf-V) 

or, more explicitly, 

'sh(C ? (±) + st7)\ sh(2# } + n) sh(5 ? (±) T a.-|) ch(f; (±) =f ei/3_ - f) 



k sh(t>f } - s?7) / sh(2t>f } - rj) sh(f;f } ± a_ + f ) ch(tff ) ± ei/3_ + £ 
x ■ 



sh(^ ±j =F 6 2 q + - f ) ch(^ =F 6 3 /3 + - g 

2> 



sh(i)f ; ± e 2 a+ + |) ch(vf 1 ± e 3/ 3 



TT sh ^f } ~ + g) sh (^ (±) + ^ + g) • , w(±) ,o 9n 

11 ,/~(±) T(±) : r /~(±) , ~(±) :> j- 1 '---'^ » 



* sh(-uj ±) - - 77) sh(t>] ±) + - 77) 



where m = Vj + rj/2. 

We have investigated the completeness of this solution for small values of s and TV 
numerically using a method developed by McCoy and his collaborators (see, e.g., [38J), 
and further explained in [21]. We find that, for k > 2sN+l (and therefore < — 1), 
all the eigenvalues of t^ ,s \u) are given by A(2>*)(-)( M ). Similarly, for k < -(2sN + 1) 
(and therefore < —1), all the eigenvalues are given by A^' s ^ + \u). Moreover, for 
\k\ < 2sN — 1 (and therefore both and are nonnegative), both A^5' s )(~)(m) 

and A^5' s )(+)(u) are needed to obtain all the eigenvalues. 

Some sample results are summarized in Tables [1]|H For example, let us consider 
Table (TJ which is for the case N — 2, s — 1. According to (I3.16p . k must be odd for 
this case. For such values of k, the number of eigenvalues of t^' s \u) which we found 
that are given by A(i s )(-)(u) and A^< S )M (u) in (13.181) are listed in the second and 
third columns, respectively. Notice that, for each row of the table, the sum of these 
two entries is 9, which coincides with the total number of eigenvalues (3 2 ). A similar 
result can be readily seen in the other tables. § These results strongly support the 
conjecture that A^' s ^~^(u) and A^' s ^ + \u) in (I3.18P together give the complete set of 
(2s + 1)^ eigenvalues of the transfer matrix t^' s \u) for \k\ < 2sN — 1. 

It may be worth noting that for \k\ = 2sN — 1 (and, therefore, either or 
M^-) vanishes), the Bethe Ansatz with vanishing M still gives 1 eigenvalue. (See again 
Tables HJH1) In the algebraic Bethe Ansatz approach, the corresponding eigenvector 
would presumably be a "bare" reference state. 

We observe that for s = 1/2, our solution coincides with the one in Our solution 
is similar to the one given by Doikou [32| 133] . but with a more general constraint 
on the boundary parameters, and with a second set of BAEs which is necessary for 
completeness. 



4 Hamiltonian for s = 1 

In this section we give an explicit expression for the spin-1 Hamiltonian, its relation to 
the transfer matrix, and its eigenvalues in terms of Bethe roots. In order to construct 
the integrable Hamiltonian for the case s — 1, we need the transfer matrix with spin-1 
in both auxiliary and quantum spaces, i.e., ^'^(u). Since we seek an explicit expression 
for the Hamiltonian in terms of spin-1 generators of su{2) (which are 3x3 matrices), 

4 For the conventional situation that the boundary parameters (a±, f3±, 9±) are finite and indepen- 
dent of N, the constraint (|3.14[) requires that also k be finite and independent of N, in which case 
\k\ < 2sN - 1 for N -> oo. 

5 It would be interesting to find a formula that would generate the entries in these tables for general 
values of N, s and k. 
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we now perform similarity transformations on the fused R and K matrices which bring 
them to row-reduced form, and then remove all null rows and columns. That is, 

Rreduced( U ) = An,02^6i,&2 -^{ai,a 2 },{6i,6 2 } ( U ) Au,a2A>i,&2 ' 



K7}dlced( u ) — An, as K^J 2 (u) A a * a2 , (4.1) 



where 



A 



( 1 \ 

\ \ 

\ -\ 

V 01/ 



(4.2) 



The reduced R and _?T matrices are 9x9 (instead of 16 x 16) and 3x3 (instead of 
4x4) matrices, respectively. It is convenient to make a further similarity ("gauge") 
transformation which brings these matrices to a more symmetric form, 

*2£L2(«) = (b ® b) R^Liiu) (s- 1 ® s- 1 ) 

*2hL£(«) = B^cedN^" 1 , (4-3) 

where 5 is the 3 x 3 diagonal matrix 

S = diag(l,-V2ch77,l). (4.4) 

We define t^ 1 ' 1 ' 9 *(u) to be the transfer matrix constructed with these R and K matrices, 
i.e., 

+(i,l) gtr \ _ t K +W gt / \ T (i,i) a* / \ K -{l) gt / a £(1,1) at / a ^ 

fc l"J — u ^reduced 0\ U ) 1 reduced 0\ U ) ^reduced 0\ U ) 1 reduced 0\ U ) ' ^- d J 

where the monodromy matrices are constructed as usual from the -Reduced ( u )' s ' anc ^ 
the auxiliary space is now denoted by "0". Since ^'^(u) and t^ 1 ' 1 ' 9t (u) are related by 
a similarity transformation, they have the same eigenvalues. Finally, it is again more 
convenient to work with a rescaled transfer matrix, 

?u , S M2 a ) s t(2„ + 2,) i(u , 

[sh w sh(u + i])\ zlv 

which, in particular, does not vanish at u — 0. 

As noted by Sklyanin [8J, the Hamiltonian is proportional to the first derivative 
of the transfer matrix 

H = cSi^v\u) +c l, (4.7) 
au «=o 
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where 

d = -chr/|l6[sh2r/shr/] 27V sh3r/sh(a_ - |) sh(a_ + |) ch(/3_ - |) ch(/3_ + |) 

x s h( a+ -|)sh(a + + |)ch(/3 + -|)ch(/3 + + |)} _1 . (4.8) 

We choose C\ so that the bulk terms of the Hamiltonian have a conventional normal- 
ization (see (14.101) ). and we choose Co so that there is no additive constant term in the 
expression (14.91) for the Hamiltonian, 

N-l 

H=Y,H n , n+ i + H b . (4.9) 

71=1 

The bulk terms H n>n +i are given by [1] 

H n , n+1 = a n - {a n f + 2 sh 2 rj [< + {S z n f + (S z n+1 ) 2 - (<) 2 ] 



- 4sh 2 (i)(^< + <<) , (4.10) 



where 

Cn = S n ■ S n+ i , = S^S^+i + S^S^ +1 , u z n = S*S^ + i , (4-11) 

and S are the standard spin-1 generators of su{2). The boundary terms H b have the 
form [15] 



# 6 = a^S*) 2 + ^S* + a 3 (Sf) 2 + a^) 2 + a 5 S? S* + a^S* 

+ a 7 S? 5+ + agS'f S? + (aj <-> and 1 JV) , (4.12) 

where iS 1 * 1 = S x ± 2iS y . The coefficients {aj} of the boundary terms at site 1 are given 
in terms of the boundary parameters (a_, /?_,#_) by 

ai = -a (ch2a_ — ch2/3_ + chr/) sh2r/sh?7 , 

a2 = -ao sh 2a_ sh 2/5_ sh 2r] , 

a 3 = — -aoe 26 - sh2?7sh?7 , 
8 

a 4 = — a e~ 26 ~ sh 2?7 sh rj , 
8 

/ 77 77 \ 3 

a 5 = aoe 6 * - ^ch/3_ sh a_ ch - + ch a_ sh/3_ sh -J sh?7 cha , 

/ 7^ Tj \ 3 

a,Q = aoe -6- ^ch/3_ sha_ ch - + cha_ sh/3_ sh -J sh?7ch2 r/, 
a 7 = —a e e ~ ych. /?_ sh a_ ch — — ch a_ sh /3_ sh — J sh 77 ch 5 77 , 

/ Tj Tj \ 3 

— aoe -6 * - f ch /3_ sh a_ ch ch a_ sh /3_ sh - ) sh 77 ch 2 " 77 , (4-13) 

V 2 2 / 
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where 



a 



Bh(a_ - \) sh(a_ + |) ch(/L - |) ch(/3_ + |) 



Moreover, the coefficients {6.;} of the boundary terms at site iV are given in terms of 
the boundary parameters (a + ,P + , 6+) by 



h = ^ 



(4.15) 



We now proceed to find an expression for the energies in terms of the Bethe roots. 
It follows from (14. 7p that the eigenvalues of the Hamiltonian are given by 

E = c 1 -^A {1 ' 1) 9t (u) +c . (4.16) 
du 



u=Q 



Furthermore, 



A (U> "(«) = ffi^AM(,), (4.17) 
|sn u sh(u + r])\ zn 

where we have used (14. 6 p and the fact that A^ 1 ' 1 ' 9t (u) = A^ 1 ' 1 ^^). ^From the fusion 
hierarchy (I2.16P with j — s — 1, we obtain (after performing the shift u — > u + r)/2) 
the following relation between AP-'^^u) and A^' 1 )^): 

A (1,1) (u) = A^\u - |) A^u + |) - + |) . (4.18) 

Recalling ([2~Tgj) 

A(3.i)( u ) =^(i.i)( w )^A(3. 1 )( u ), (4.19) 

and also our result (13.181) for A^' 1 )(u) ) we finally arrive at an expression for the energies 
in terms of the Bethe roots 

E = sh\2 v ) V m „ + N - s] + C W , (4.20) 

^ sh(^ (±) - rj) sh(^ (±) + n) V shr/ ^ 

where , uj ± ' ) = + rj/2 as in (13.211) , and are constants whose cumbersome expres- 
sions we refrain from presenting here. (These constants are independent of N, but do 
depend on the bulk and boundary parameters and on {e^}.) 

We have verified that the energies given by the Bethe Ansatz (I4.20p coincide with 
those obtained by direct diagonalization of the Hamiltonian (14 .9p for values of N up to 
4. For the case N = 3, some sample results are presented in Tables [5], El for boundary 
parameter values corresponding to k — 1. Hence, as already noted in Table [2} 17 levels 
are obtained with A^> 1 )(~)(m) and are listed in Table EJ and 10 levels are obtained with 
A ( i 1)(+) (u) and are listed in Table [61 Together, they give all 27 energies obtained by 
direct diagonalization of the Hamiltonian. 

For s > 1, it is in principle possible to proceed in a similar way. However, the 
computations become significantly more cumbersome, and we shall not pursue them 
further. 



11 



5 Discussion 



We have found a Bethe Ansatz solution of the open spin-s XXZ chain with general 
integrable boundary terms (13 . 1 8[) - (13.211) , which is valid for generic values of the bulk 
anisotropy parameter t], provided that the boundary parameters satisfy the constraints 
(13.141) . ( 13.161) and ( 13.171) . We have presented numerical evidence that this solution is 
complete, and we have explicitly exhibited the Hamiltonian and its relation to the 
transfer matrix for the case s = 1. The apparent correctness of this solution provides 
further support for the validity of the identification (13.11) of the Q operator as a transfer 
matrix whose auxiliary-space spin tends to infinity. 

A drawback of this solution is that for \k\ < 2sN — 1, one does not know a priori 
in which of the two "sectors" (i.e., A^' s ^~)(u) or A^< s )(+)(u)) a given level - such 
as the ground state - will be. For the case s = 1/2, alternative Bethe Ansatz-type 
solutions have been found [31] which do not suffer from this difficulty, and for which 
the boundary parameters do not need to obey the constraints (13.141) , (13. 16[) and (I3.17P . 
However, those solutions hold only for values of bulk anisotropy corresponding to roots 
of unity. (For yet another approach to this problem, see [28].) Perhaps such solutions 
can also be generalized to higher values of s. 

Part of our motivation for considering this problem comes from the relation of 
the s = 1 case to the supersymmetric sine-Gordon (SSG) model [39], in particular, 
its boundary version [40]. Indeed, Bethe Ansatz solutions of the spin-1 XXZ chain 
have been used to derive nonlinear integral equations (NLIEs) for the SSG model on 
a circle [JT] and on an interval with Dirichlet boundary conditions [12]. With our 
solution in hand, one can now try to derive an NLIE for the SSG model on an interval 
with general integrable boundary conditions, and try to make contact with previously 
proposed boundary actions and boundary S matrices |40j . 
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A Conjectured relations for the fusion hierarchy 



One way to derive the fusion hierarchy (12.161) relies on the existence of the relations 

2j-i 3 

-^12...2j-l -^12...2j -^12...2j-l = ^12...2j-2^2j-l,2j + ^k,k+l ^ ^ ^k,k+l ' J = 2' ^' ' ' •C^'l) 

fe=l 

where P\ 2 ___ n = 1 — -PiL.m f° r some 

set of matrices {X^ , X^~^}. To our knowl- 
edge, these relations have not been proved for general values of j. However, we have 
verified them up to j = 3. In particular, for j = 3/2, 

X (1) = I (v 23 - V l3 ) , X« = l - (v 12 - V V3 ) . (A.2) 

This result is equivalent to an identity found by Kulish and Sklyanin, see Eq. (4.15) 
in the second reference of [2]. For j = 2, we find 

X {1) = ^{V 23 -V 13 + V 2A -V 1A + 2V l3 V 2A -2V 14 V 23 ) , 
X {2) = -(P13-P12 + P34-P24 + 2P12P34-2P13P24), 

X {3) = \{V 13 -V u + V 23 -V 2i ). (A.3) 
Moreover, for j = 5/2, 

X (1) = |p 15 p 14 + p 15 p 24 - ^p 15 p 34 + Ip 35 p 14 
5 5 2 

- ^35^13 " ^45^14 - ^45^13 + ^V M V 13 , 

10 10 o 15 

X {2) = \v ib V 2A + ^-V 15 V 13 - \v 2h V U + ^25^24 + 4^25^34 

5 15 5 5 15 

- ^V 25 V 13 - ^V 45 V 2A - ^-V i5 V 13 + ^-V M V 13 , (A.4) 
15 15 15 15 

x {3) = —V l5 V u -—V 15 V u + —V 25 V li +—V 25 V M 

- ^v 35 Vu - \v 3b v 3A - -v 35 v 12 , 

x {4) = ~v 15 v u + lv 15 v 2A + lv 15 v u -lv 15 v 23 . 

5 5 5 5 

We omit our lengthy results for j = 3. We emphasize that the matrices {X^} are by 
no means unique: for instance, only the antisymmetric part X[ k k+1 2 ■ — x[ k+1 k 2 - 
matters in the expression of X^ k \ In fact, for the purpose of deriving the fusion 
hierarchy, the explicit matrices are not important, as the terms ^ k P k + k+1 X^ k > P k ~ k+1 
in (lA.ip do not contribute. What is important is only the fact that such matrices exist. 
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We also observe that the symmetrizer P+ n (12. 3p can be expressed as a sum of 
products of commuting permutations: 

[n/2] 

A + ..n = + E °* E CT ( 2 ) V °$) CT ( 4 ) ' ' ' ^(2<-l) »(20 ( A - 5 ) 

where <S n is the permutation group of n indices. The coefficients are given by the 
recursion relations 

4 n+1) = ^(al?>-n\aM) (A.6) 

>+ x ) - 1 f 1 „(») i /) „(») o„(n) 



(n+1) 31 + 1 (n) _ l + l (n) n + 2 - 21 (n) 

^ " (n + iy e ~(n+l) 2 + (n + lf ^ ' " " 



^1 

- (A.i 
2J v 



(n+1) _ n-2[n/2] (n) 

V/2]+l - (n + 1) 2 ^ 1 ' 

with a|, = 1 and the convention af 1 ^ = when £ > [n/2]. Above, [ ] denotes the 
integer part. Note that a!jp = 1 is consistent with our previous convention Vn = 1, 
which implies P x + = 1 and Pf = 0. From the recursion, one then finds e.g. ag = | and 
= |, which reproduces the result P^ = |(1 + P^)- Let us stress that, because of 
the sum on all permutations in S n , a term containing I permutations has a multiplicity 
l\ [n — 2£)\ 2 £ with respect to a 'reduced' expression, where all the terms appear just 
once. Identity ( 1A.5I) may help in the proof of the conjecture (1A.1I) . 

One can prove relation ( 1A.5I) by recursion. Starting from the relation at level n, 
the recursion relation 

(n + 1) Pt.n + 1 = \ 1 + E ) P l + -.n = P L..n ( 1 + E ^ 



£=1 



and the identity (valid for any spaces a, 6, c) 

Pa6 fPa6 + Vac + ?0 = V ab + P ac + P 6c 



show that the relation is also obeyed at level n + 1. A careful analysis of the different 
terms leads to the relations (IA.6HA.9l) . 
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k 


no. given by A^^X \u) 


no. given by A^2' 1 )(+)(it) 


5 


9 





3 


8 


1 


1 


6 


3 


-1 


3 


6 


-3 


1 


8 


-5 





9 



Table 1: Number of eigenvalues of t^' s ^(u) given by A^ s )(±)( n ) for N = 2, 
s — 1. The total number of eigenvalues is (2s + 1)^ = 3 2 = 9. 



k 


no. given by A^^X )(u) 


no. given by A^2' 1 K+)( lt ) 


7 


27 





5 


26 


1 


3 


23 


4 


1 


17 


10 


-1 


10 


17 


-3 


4 


23 


-5 


1 


26 


-7 





27 



Table 2: Number of eigenvalues of td ,s \u) given by A(2< S )(±)( M ) for N = 3, 
s = 1. The total number of eigenvalues is (2s + 1)^ = 3 3 = 27. 
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k 


no. given by A^'a^ ^(u) 


no. given by A^'^Wfu) 


7 


16 





5 


15 


1 


3 


13 


3 


1 


10 


6 


-1 


6 


10 


-3 


3 


13 


-5 


1 


15 


-7 
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Table 3: Number of eigenvalues of t^ ,s \u) given by A(2> S )(±)( M ) for N 
s = 3/2. The total number of eigenvalues is (2s + 1)^ = 4 2 = 1 



k 


no. given by A^'i^ \u) 


no. given by A^'^+^it) 


10 


64 





8 


63 


1 


6 


60 


4 


4 


54 


10 


2 


44 


20 





32 


32 


-2 


20 


44 


-4 


10 


54 


-6 


4 


60 


-8 


1 


63 


-10 





64 



Table 4: Number of eigenvalues of t^ ,s \u) given by A ( 5' S )( ± )( M ) for iV 
s = 3/2. The total number of eigenvalues is (2s + 1)^ = 4 3 = 6 
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E 


Bethe roots Vj 


-8.78796 


0.0781924 ± 0.150582 i , 0.573709 


-7.99601 


0.377364+1.5708 i, 0.0718753 ± 0.150316 i 


-5.5443 


0.191917 ± 0.145165 i , 0.529223 


-5.07143 


0.375505+1.5708 i , 0.164075 ± 0.148506 i 


-4.31229 


0.158193 ± 0.299905 i, 0.158448 


-3.36195 


0.166008, 0.717455 ± 0.259354 i 


-2.87198 


0.358903+1.5708 i, 0.156172 , 0.784233 


-2.86245 


0.33466 ± 0.286332 i , 0.337633 


-2.69332 


0.371101+1.5708 i , 0.292713 ± 0.157296 i 


-2.31742 


0.290731 + 1.5708 i, 0.617492+1.5708 i , 0.146609 


-1.52379 


0.484424, 0.621449 ± 0.318594 i 


-1.18428 


0.356639+1.5708i, 0.464322, 0.659312 


-0.780678 


0.288176+1.5708 i, 0.610874+1.5708 i , 0.368261 


-0.379026 


0.879352 ± 0.483137 i , 0.944398 


-0.0249248 


0.337585+1.5708 i, 0.934391 ± 0.266345 i 


0.389221 


0.277491+1. 5708i, 0.580415+1.5708 i, 0.97389 


0.838091 


0.245314+1.5708 i, 0.477481+1.5708 i, 0.814847+1.5708 i 



Table 5: The 17 energies and corresponding Bethe roots given by A^ ,1 )( 

ioi N = 3,s = l,k = l,r) = 0.3i,a_ = 0.7i,/3_ = 0.2, 6_ = 
0.5i,a+ = l.2i,(3+ = -0.2,6+ = — l.li , {e^} = +1. 



E 


Bethe roots 


-9.55066 


0.0900396 ± 0.151265 i 


-5.71507 


0.244797 ± 0.132886 i 


-4.09573 


1.50013 i, 0.182899 


-3.74447 


0.786256 i, 0.174481 


-3.02558 


0.514399 ± 0.220939 i 


-2.07231 


1.48515 i, 0.620007 


-1.79241 


0.80571 i, 0.568604 


-1.1462 


0.350408 +1.5708 i , 1.38463 i 


-0.791024 


0.093628+1.5708 i, 0.842728 i 


-0.634944 


0.979155 i, 0.863041 i 



Table 6: The 10 energies and corresponding Bethe roots given by h.^ ,l ^ + \u) 
for N = 3,s = l,k = l,r] = 0.3z,a_ = 0.7z,/3_ = 0.2, 6_ = 
0.5i,a+ = 1.2i,(3+ = -0.2,6+ = -l.li , {ej = +1. 
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